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We evaluate the distribution function of mesoscopic fluctuations of orbital magnetic response in 
finite-size level-quantized metal particles and Aharonov-Bohm rings for temperatures smaller than 
the mean level spacing. We find a broad distribution with the reduced moments much larger than 
the mean. For strong spin-orbit interaction we find very long tails due to thermal activation of large 
effective moments of the electrons at the Fermi level. 
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I. INTRODUCTION 

The question of mesoscopic fluctuations of orbital magnetic response has been extensively studied within the 
. perturbation theory approach which, generally, holds when the temperature and/or level broadening is larger than 
r/2 ' the mean level spacing A jij . In 2D (which we only consider here) , the variance of the magnetic moment was predicted 
<D ' to be"' 



<MO ~ A (k F ef (£) 2 in (J) ~ A (k F ef (^) 2 In (J) (1) 



where 



T3 . 

T*=max{T,T if 1 } (2) 

O ■ 

with the following notations: /is is Bohr magneton, kp is the magnitude of the Fermi wave vector, I is the electron 
mean-free-path, </> is the flux through the sample or the Aharonov-Bohm (AB) flux for rings, </>o = 27r/e is the flux 
quantum in units where U = c = 1, E c ~ D/L 2 , where D = vf£/2 is the diffusion coefficient, and L the sample size 
^ . (typically a ring or disk circumference), and 



7 N < I {^-) - ^^(k,.n (8a) 



The logarithmic dependence changes its form when 



rp \ 1/2 



(4) 



The significance of the flux scale 4> c is such that at T ~ A it determines the scale of transition from Gaussian Orthogo- 
\ nal (GOE) to Gaussian Unitary (GUE) Ensemble. It is also the scale of linear response in the problem. (While eq. (jl]) 
indicates that, aside from logarithmic corrections, the linear response regime extends to <f> ~ c/)q for the fluctuations, 
it is obtained in the so called "zero-mode" approximation [Q for systems with simply connected geometries. However, 
the higher modes may need to be summed up to correctly evaluate the fluctuations [Q . Consequently, unless specifi- 
cally stated otherwise, we shall limit our consideration to <f> < <j) c where the response is linear both for the fluctuations 
• i-h . and the mean, for simply connected and AB geometries.) 

£ : 

_ _ i 

1 In an earlier work, Altshuler and Spivak considered current fluctuations in SNS junctions with broken time-reversal 
symmetry. This problem is equivalent to the persistent current fluctuations in AB rings. If the time-reversal symmetry is 
not broken and the Cooperon contribution is taken into account, their result would be equivalent to Ref. jjj], including the 
logarithmic term. 

2 Insofar as the exact numerical coefficient in eq. (|f|), both Refs. [[jj and j| used the Euler-Maclaurin method for Matsubara 
summation with a resulting integral over the continuous variable x = mT/E c , dx — T/E c <C I. However, for small x a more 
accurate procedure must be used. It is explained in Ref. || for the mean response (Appendix II) and can be trivially extended 
to fluctuations. 
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The average magnetic moment was predicted to be 



W~ (5) 



which, for T ~ A, becomes 



^ (k F £) £ ~ n B ejE (fat) , $ < <f, c 

Comparing (|J) and (|l|) on the scale < C , we see that (Af) and (JA/ 2 ) 1 ^ 2 are of the same order of magnitude. While 
numerically the latter is larger Q , it is clear that the diamagnetic response is much less likely than the paramagnetic 
response. 

In this work we address the regime in which level-quantization becomes important, namely, T <C A. In this case, the 
perturbation theory approach, which uses the level density correlation function and standard thermodynamics, is no 
longer applicable and one needs to use a single electron picture in conjunction with thermal occupancy of a two-level 
system. In previous papers, we considered the mean response of GOE and Gaussian Symplectic Ensemble (GSE), 
the latter being the case for strong spin-orbit (SO) interaction. We have argued the for GOE the mean response can 
be formulated in terms of a single-electron van Vleck response at the Fermi level || and for GSE in terms of the 
effective magnetic moments/persistent currents of electrons in the last occupied (Fermi) state Here, we apply 
these approaches to mesoscopic fluctuations. 



II. GOE 



As explained in Ref. || , the mean orbital response can be understood in terms of the van Vleck (vV) response that 
involves virtual transitions from the last occupied level to the first unoccupied level Sf of the Fermi sea 



[i\M z \f) 



H 2 



if 



H 2 



£ f 



(J) 



where M z is the magnetic moment along the magnetic field H (perpendicular to the sample) . There are three possible 
sources for non-self-averaging (fluctuations) based on this picture. The first two, based on eq. @), are the fluctuations 



of 



M, 



if 



and the distribution of the level spacing at the Fermi level (that is, the distribution of the values in the 

denominator of (Q)). We expect the latter to be dominant and will neglect the former. The third, which may lead to 
occasional diamagnetism, is the system-dependant nature of the Fermi-sea cancellation between the diamagnetic and 
paramagnetic contributions to the total magnetic response 5), |m[; it is not studied here. 



The mean value of the vV response is obtained as follows. 



M, 



if 



can be found using the scmiclassical approximation 



(and, more precisely, using the result for the magnetic dipole absorption) || and is given by 



A l, 



if 



In what follows, we will use a dimensionless measure x for the level spacing at the Fermi level, 

Sf -Ei = xA 

Averaging with the GOE Q distribution function for the nearest energy levels, we find |] 



M, 



if 



H —— 

2A 



\ XL \H 2 A(k F £) 



exp 



irx 



7TV 

ax = 

2 



M, 



if 



H' 2 



(8) 
(9) 

(10) 
(11) 



Here s is the level degeneracy (s = 2, on the account of spin), v — sA _1 is the density of levels, xl is the Landau 
susceptibility and A is the sample area. As was pointed out earlier, it is assumed that r^ 1 <C A, with the opposite 
limit corresponding to GUE. 
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To consider the fluctuations, we notice that the second order perturbation theory is used in derivation of (Q). 
Consequently, the limit of its applicability jlo| is, per ((?]) and 



€vV 



x 



< xA 



that is 



x > r 



s x /a) 



1/2 



e/ - e 2 » (r H 



1/2 



(12) 



(13) 



An important observation about the above conditions is that the range of linear response is sample-specific and depends 
on the Fermi level spacing of a particular particle. On the other hand, when describing the statistical distribution of 
the magnetic energies, it is only meaningful to consider their values at a given field (or a range of fields). 
An estimate of the fluctuation can be obtained by substituting the r.h.s. of ( fl3| ) into eq. ((T 



1/2 



(14) 



Below we will show that (M])) corresponds to the result for the reduced higher cumulants for the distribution of 
magnetic energies and is due to the systems with small level spacing at the Fermi level. 

Since the reduced moments grow with the order of the moment, the cumulants and the moments should be of the 
same order of magnitude and it is sufficient to evaluate the latter. We begin with the evaluation of the variance/second 
moment of the vV response 



h 



dx 

— exp 



TTX 



1 H 



In 



A 



(15) 



(omitting subscript "vV"). The higher moments and the reduced moments can be evaluated similarly and are given 
by 



H 



-n+1 



^)" 2 (^- 



1/2 



TTX 



1/n 



dxexp [ — 



(16) 
(17) 



in agreement with the estimate (14). Alternatively, eqs. (|15|)-(|17|) can be evaluated using the distribution function 
for the fluctuations which is found as 



P(e) 



dxS e 



JL 

x 



x exp 



TTX 



'I 'H 

2 e 3 



exp 



' H 



4e 2 



e < (r^A 



1/2 



(18) 



with the case e > (r if 1 A) 1 2 (or, equivalently, x < (r i/ 1 /A) 1 2 ) addressed below. 

A comparison should be made between (|l5|) and the perturbative result ([!]). First, assuming for the latter that 
T ~ A (the limit of applicability of the perturbative approximation) , we find that the corresponding energy fluctuation 



tpert 



(19) 



The difference in log factors is because (|l5| ) is obtained in a two-level approximation, whereas (19) takes into account 
contribution of levels within E c of the Fermi level; otherwise, the results are in qualitative agreement. (The mean 
values e v y are also in agreement A). Although the higher moments had not been evaluated perturbatively, it is 

expected that n \/£p er t ~ an d are smaller than obtained in a two level picture. This should be anticipated since 



the two-level response is very sensitive to the variations of the energy level spacing at the Fermi level. 

1/2 

For x < (r^ 1 /A) the perturbation theory evaluation of e v y (0) is no longer valid. This applies to the fraction 
~ r if 1 /A of all systems with sufficiently small Fermi level spacing. There are two possible approaches to this case 
that yield, essentially, the same result. In the first approach, the Fermi level and the first unoccupied state can be 
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viewed as, effectively, a doubly degenerate state with respect to the magnetic field perturbation. Using the secular 
equation 10], we find that the magnetic field will split the levels in the amount 



M, 



if 



V~(r^A) 1/2 (20) 



(The second order term is still given by eq. (m) but with the final state / being the next unoccupied state at a distance 

1 /2 — 

greater than (r^ 1 A) from the Fermi level.). Alternatively, one can argue that such systems are effectively in the 
GUE regime. The perturbative expression for the mean value is 



1 /2 

(pert (GUE) ~ (t# A) (21) 



1 /2 

which is obtained from the second eq. using the substitution A — > (r i? 1 A) , where the latter is the upper 
bound of the level spacing in the systems that have effectively undergone the GOE — > GU E transition. 

The physical interpretation of the above results may be as follows. In a particular systems, the linear response 
regime is determined by the Fermi level spacing in this system, as per eq. (|l2|), 

M~glgg*3 (22) 
xA v ; 

As a result of linear response, the magnetic moment approaches the value 

M ~ /is (fc^) 1/2 (23) 

with little change for larger fields (until r^ 1 approaches A - see below). The statistical distribution of the response 
implies the need to consider the same range of fields for all systems. The mean magnetic moment is given by 

{M) ^fM^l^% B{kF ^ (24 ) 

The slowly decaying distribution function 

P( £ )(x^,r H 1 < e <(r H 1 A) 1/2 (25) 

is due to systems with small Fermi level spacing, where the magnetic moment may become as large as (p3|). 

At Tfj 1 > A, the single-level ansatz is no longer valid and the results of Ref. (as described in the Introduction) 
should be applied instead; at t^ 1 ~ A the latter are consistent with the present results. 



III. GSE 



It was argued |Tl]], H that in a GSE the electrons at the Fermi level have a magnetic moment 

M ~ [Ib (k F £) 1/2 (26) 



(which is, incidentally, of the same order of magnitude as given by (|23j)) induced by the spin-orbit interactions. For 
a system with even number of electrons, the magnetic moments of the two Fermi-level electrons cancel each other. A 
system with an odd number of electrons, on the other hand, should have a permanent magnetic moment || and a 
Curie- like susceptibility. This is in incomplete analogy with the purely spin magnetism Jl2| , the sole difference being 
that M — /is in the latter case. Consequently, we refer to Ref. Q for the discussion of the distribution function of the 
susceptibility obtained in the two- and three-level approximation for the even- and odd-electron systems respectively; 
in both cases, one finds long tails due to thermal activation of the magnetic moments. 

It should be noted that the single-electron magnetic moment (|2^) can be attributed to the large electron g-factor, 

g ~ (kpi) 1 ^ 2 , induced by the spin-orbit interaction 0, [Of, pi]. Furthermore, it was argued |13|, p| that the latter 



can fluctuate by as much as the order of magnitude. We point out, however, that in the temperature regime T <C A 
considered here, the thermal activation effects should be dominant as they may produce much larger fluctuations. 
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IV. DISCUSSION 



The two-level picture considered here predicts the large fluctuations of the orbital magnetic response due to the 
possibility of small level separations at the Fermi level. Obviously, the possible observation of these effects imposes 
restrictions on the temperature and the magnetic held. Moreover, such restrictions are sample specific since the 
magnetic energy and the temperature must be compared with a particular Fermi level separation of the nearest 
states. (See, for instance, eqs. ( |l2|) and (|l3|)). On the other hand, a statistical description of the fluctuations implies 
a common (range of) magnetic field and temperature for all systems. Consequently, while in a given system one can 
observe the linear response regime due to the vV energy (Q) by sufficiently reducing the magnetic field, the Fermi-level 
spacing in this system may be smaller than the magnetic energy corresponding to the field at which all systems are 
analyzed. The latter is critical for understanding of the effective cut-offs in the analysis of the distribution function. 

2 

and the details of 



if 



The issues that remain to be understood are the fluctuations of the matrix elements 

the cancellation between the van Vleck paramagnetism and the precession diamagnetism over the Fermi sea. These 
will require further, largely numerical, studies. Also, our formalism should be applicable, with minor changes, to 
orbital magnetism of integrable systems, such as a rectangle with incommensurate sides Jl5[ . We hope to address this 
problem in a future work. 
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